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RATIO ; Detailed Current Edit History 
POWJJ = longword to power longword giving longword result 


N 5 
OTS$POWJJ - INTEGER®4 ** INTEGER*®4 power routine 16-SEP-1984 AX/VMS Mac £8 v04-00 Pa 1 
1-005 , -SEP=-1984 4: 38 $§ YRTHRTL. Sh SRCJOTSPOWJJ MAR; 1 ad (1) 
TITLE QISSA QUI = INTEGER*4 ** INTEGER*4 power routin 
~IDENT /1-005/ 3; File OTSPOWJJ.MAR Edit: SBL1005 


MPARABRARASRASALALAALALASAAAAELESASASER REARS EERE ARERR REESE SERRE REESE SEER REESE SS OS 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGIT AL EQUIPMENT C0 CORPORATION, MAYNARD, MASSACHUSETTS. 


RIG 
TH fo’ ‘Keche iY FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
0 IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
SION OF THE ABOVE Slt had NOTICE. THIS SOFTWARE OR ANY OTHER 
S_THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 
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** ® 
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** © 
‘# QO * 
;* TRANSFERRED. © 
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** 
°* * 
** " 
** © 
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** 7 
°* ra 
** * 
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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
oeePORATI Ow. NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONS 
SOFTWARE ON EQUIPMENT WHIC 


ITY FOR THE USE OR _ RELIABILITY OF ITS 


BILIT 
IS NOT SUPPLIED BY DIGITAL. 
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; FACILITY: Language support Library - user callable 
ABSTRACT: 
Integer longword base to integer Longword power. 


Integer overflow can occur if the result exceeds a lLongword. 
Undefined exponentation can occur if base is 0 onl power is 0 or negative. 


VERSION: 0 
HISTORY: 
; AUTHOR: 
Thomas N. Hastings, 5-May-77: Version 0 


MODIFIED BY: SUSAN HUBBARD AZIBERT 
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ORY ; Detailed Current Edit History 6-SEP- MTHRTL.SRCJOTSPOWJJ.MAR; 1 
-SBTTL HISTORY ; Detailed Current Edit History 
4 ; Edit History for Version 01 of OTS$POWJJ 
5 ; version 04 = changed module name to forpowJJ 
+4 § 3 version § = changed error routine from PATHSERROR to MTHSSERROR 
00 3 version - removed W* on call to MTHS$S$ stlana’s seve code with MOVZBL. 
000 8 ; gers ton 8 = change MTHSSERROR to gps ag ee! 
000 9; O-11. = fix case instruction 8, -Feb-7 
00 60 ; 1-001 - Y ate version number and coprr ht ry ges ne NOV-78 
009 61: 1-008 = Change, ATH. UNDEXP to MTHSK UNDEXP. rp 
4 6 3 1- 3 - o the J girective. of S 
0 65 ; 1-004 - dec late MB yo May-19 
0 64 ; 1-005 = Use general mode Seerenbine. SBL , etapa 
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«SBTTL DECLARATIONS 


> INCLUDE FILES: 


> EXTERNAL SYMBOLS: 


-DOSABL GBL 
~EXTRN MTHSK_UNDENE 
-EXTRN MTHSSSTGNAL 


: MACROS: 


> EQUATED SYMBOLS: 
base = 4 
exp = 8 


OWN STORAGE: 


: PSECT DECLARATIONS: 


ODWNAVE WN HO ODNOA UE 0 ODOUR OOOO 


OOOO 0000000000000 0009 090909 NI NIN NINN SN NIOOOAO 


16-SEP-19 
e-SEP-19 


Be 14:28:39 


; Math error routine 


> base input formal - by-value 
; exponent intpu formal - by-value 


-PSECT _OTSSCODE PIC,SHR,LONG,EXE ,NOWRT 
3; program section for OTS$ code 
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-SBTTL OTSS$POWJJ = Longword to power longword giving longword result 


;+¢ 

3; FUNCTIONAL DESCRIPTION: 

Signed longword result = Signed longword base ** signed longword exponent 
s 


® 
OW lo 
' 
1 
108 
106 ; 
: ; The signed longword result given by: 
103 ; base exponent result 
14 ; any >0 product (base * 2**i) where i is each | 
8 16 : non-zero bit position in exponent 
0 114 ; > s 1 
$0 115 ; = = 6 Undefined exponentation 
9000 116 <0 1 | 
38 118 ; >1  «€0 0 
00 119; = 1 < 0 1 
000 120; = 0 <0 Undefined exponentation 
44 121 ; = -1 < 0 and even 1 
000 «(1 § 3 = =] < 0 and odd =1 
Bs te bo pee | 
0000 125; Integer overflow can occur. 
0000 =—1 8 5 Undefined exponentiation occurs if base is 0 and 
0000 127 ; exponent is 0 or negative. 
0000 128; 
0000 129 ; CALLING SEQUENCE: 
0000 130; 
444 } 1 ; Power.ww.v = OTSS$POWJJ (base.rw.v, exponent.rw.v) 
0900 5 > INPUT PARAMETERS: 
0000 135 ; 
444 : $ ; states arti 
000 138: 
$38 133 : Serres Seen meet 
900 141: | 
44 196 : SRPLECEY QUTPUTS: 
9 144 ; 
102 : FUNCTION VALUE: | 
0 9 147 : Longword integer base ** exponent 
0 145 > SIDE EFFECTS: | 
134 ; SIGNALS SS$_ARITH with integer overflow hardware code if | 
1 § 3 snteger overflow. 
1 3 SIGNALS MTH$_UNDEXP (82 = * UNDEFINED EXPONENTATION') if 
: 3 : base is 0 and exponent is 0 or negative. | 
156 3-- | 
139 | 
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TSSPOWJJ = INTEGER®4 ** INTEGER®*4 power routine 16-SEP-1984 02:02:42 VAX/VMS Macro v04-00 Page 5 
#885 OTSS$POWJJ = longword to power longword g orSFb robe 99:98:45 MTHRTL.SRCJOTSPOWJJ.MAR; 1 . (4) 
13 | 
4004 160 -ENTRY OTSSPOWJJ, “M<IV, R2> i enable integer overflow 
50 D 161 MOVL » RO ; RO = initial result 
52 OB AC OD 196 MOVL exp(AP), R2 ; R2 = exponent 
1 ; 197 BLEQ EXPLEQ ; branch if exponent =< 0 
165 ;¢ 
198 3; Exponent > 0. 
167 ; Scan each exponent bit from ions. squaring base each time thru nf 
008 188 3; For each 1l=bit in exponent, multiply current base into partial result. | 
B 170° 
51 04 as D0 0008 171 MOVL base(AP), R1 ; R1 = base 
0B 5 3 OOF 7 BLBS R2, PARTI 3; branch if exponent is odd 
52 52 FF er 1g 173 SQUAR: ASHL #-1, R2, R > R2 = exponent/2 
51 1 ¢4& 001 174 SQUAR1: MULL R1, R1 ; RI = current power of base 
O1A =175 ; integer overflow will trap 
O1A 178 ; and SIGNAL SS$_ARIT 
F552. E€9 OOIA 17 BLBC R2, SQUAR ; loop if next bit in exponent is 0 
Baie 178 PARTIAL: 3 next bit in exponent is a 1 
50 51 C& O001D 179 MULL R1, RO ; RO = new partial result 
0020 180 3; integer overflow will trap 
0020 181 ; and SIGNAL SS$_ARITH | 
$2 f2 tte 7 83 2 135 ASHL #-1, R2, R2 3 = exponent/ 
FO 12 00 18 BNEQ SQUAR1 ; loop if more exponent bits are 1 | 
04 0027 184 RET : return, RO<15:0> = base ** exp 
0058 185 : RO<31:16> = 0 | 
0028 186 
88 187 ;+ 
0 188 ; Exponent is =< 0. 
00 189 ;- 
0028 190 
08 19 0028 =191 EXPLEQ: BLSS EXPLSS ; branch if exponent < 0 
002A 136 
O02A 193 ;+ 
Q02A 194 ; Exponent is = 0. 
Q02A 195 ; Undefined exponentation if base = 0 too, else return 1 
QOcA 196 ;- | 
O02A 19 
51 04 AC 00 002A 198 MOVL base(AP), R1 ; R1 = base | 
1B 13 00 199 BEQL UNDEF INED ; undefined if base = 0 too | 
18 11 00 00 BRB POWJ JX 3 return with result = | 
84 3 since base ** 0 = 1 
08 i+ 
0 04 ; exponent =< 0. 
05 ; Result is given by the following table: | 
3 Bese Result 
§ : -1 1 or -1 depending on exponent being even or odd 
0 9 : 9 yndet ines exponentation 
12: >1 0 | 
Baie 
14 
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#805 OTSSPOWJJ = longword to pover longword g 6=-SEP=-1 1984 94 $8 %$ ERTHRTL. SR SRE fOTSPOWJJ. MAR; 1 . 
15 EXPLSS: 
O2 FFFFFFFF 8F 04 AC 3 1g 108 Chet tte tm #-1, #2 5 fase. on votue f bage te ai 
: ‘ 1-1 § return = -1 or ependin 
D 3; on 8 exponent being odd or even . 
4 po aay “hs § gH veo hha exponentation 
. - § u 
i} CLRL RO ; i< =j or > #1): return RO = 0 
rh RET 
44 MINUS]: BLBC Ree POWJJX ; if exponent is even. return RO = 
MNEGL #1, RO : else return RO = -1 
POWJJX: RET ; return 
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Undefined exponentation error - 0**0 or 0*#*(-n) 
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48 
48 UNDEFINED: 
50 D4 48 CLRL RO ; if error, return result = 0 
7E 00° 8F GA 4D ir pla UNDEXP, =-(SP) ; FORTRAN eres g 

00000000'GF 01 FB 51 CAL , G*ATHS$SIGNAL 3 convert to 32-bit condition code 
28 3 ane ao hay UNDE XP 
5 3 end arg not needed since no 
2 ; "e6 *OTSSPOUR is possible. 

04 5 RET 

59 
59 
59 END 
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OTSSPOWJJ = INTEGER*4 ** INTEGER*®4 power routine y -SEP-1984 
Symbol table , Pa 198 04: 98 $$ t 
BASE 
EXP 
EXPLEQ 1 
ar 1 
MINUS 1 
MTHSSS1IGNAL ah 
MTHSK_UNDEXP eeereeee 
OTS$POwJJ 0000000 Re 1 
PARTIAL 00001D 1 

JJ 44 A 4 1 
SQUARI B000017 RT 
UNDEF INED §o080048 fe 1 

¢eeceecen ewer ere aan + 
: ; Psect synopsis : 

PSECT name Allocation PSECT No. Attributes 
OS . 000000 0 <¢ 0.) 69 
-OTSSCODE 00000059 («¢ 89.) 01 1.) PIC USR CON REL LCL 


gueceeecce 


Phase Page faults CPU Time 
Initializati ° 
Command processing 136 3 00:09 00:4 ‘ 3 
pyabol table sort 0 0:00: §°3 
Pass 2 54 0:00: 0:5 
Symbol table output 2 0:00:00.0 
Psect synopsis output 2 0: 8: 0.0 
Cross-reference output ° 0: 9-9 
Assembler run totals ey 00:00:01. 


ihe working set Limit was 750 pages. 
7 bytes (6 pages) of virtual memory were used to buf 
¢ re were 10 pages of symbol table space allocated to 


AX/VMS M 
MTHRTL. SRC 


0 V04-00 
JOTSPOWJ. MAR; 1 
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( Q.) NOP IC “USR CON ABS LCL ag woene NORD NOWRT NOVEC BYTE 
( SHR RD NOWRT NOVEC LONG 


eeeceeeoeeooocoen} 


0:00:01.48 


:00: 
700: 
:00: 
700: 
:00: 
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fer the intermediate code. 
hold 13 non-local and 1 local symbols. 


$6 source lines were read in Pass 1, producing 11 object records in Pass 2 


pages of virtual memory were used to define 0 macros. 
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Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB;2 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 


A Ho 080s 02: z | 1 
GARsit Aacro Run Statistics INTEGERS4 ss INTEGER*G power routine 19-3 eOn18Re 99:58:96 PATA ekesorsecuss.man;1 P28* By - 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:OTSPOWJJ/OBJ=OBJ$:OTSPOWJJ MSRC$:OTSPOWJJ/UPDATE=(ENHS:0TSPOWJJ) 


| 0265 AH-BT13A-SE ) NT CORPORATION 
gave VAX/VMS V4.0 D PROPRIETARY 


